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observed difference of velocities of neutrals and ions 
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ABSTRACT 

Theoretical and observational studies on the turbulence of the interstellar 
medium developed fast in the past decades. The theory of supersonic magne- 
tized turbulence, as well as the understanding of projection effects of observed 
quantities, are still in progress. In this work we explore the characterization 
of the turbulent cascade and its damping from observational spectral line pro- 
files. We address the difference of ion and neutral velocities by clarifying the 
nature of the turbulence damping in the partially ionized. We provide theoret- 
ical arguments in favor of the explanation of the larger Doppler broadening of 
lines arising from neutral species compared to ions as arising from the turbulence 
damping of ions at larger scales. Also, we compute a number of MHD numerical 
simulations for different turbulent regimes and explicit turbulent damping, and 
compare both the 3-dimensional distributions of velocity and the synthetic line 
profile distributions. From the numerical simulations, we place constraints on the 
precision with which one can measure the 3D dispersion depending on the tur- 
bulence sonic Mach number. We show that no universal correspondence between 
the 3D velocity dispersions measured in the turbulent volume and minima of the 
2D velocity dispersions available through observations exist. For instance, for 
subsonic turbulence the correspondence is poor at scales much smaller than the 
turbulence injection scale, while for supersonic turbulence the correspondence is 
poor for the scales comparable with the injection scale. We provide a physical 
explanation of the existence of such a 2D-3D correspondence and discuss the 
uncertainties in evaluating the damping scale of ions that can be obtained from 
observations. However, we show that the statistics of velocity dispersion from 



^Nucleo de Astrofisica Teorica, Universidade Cruzeiro do Sul - Rua Galvao Bueno 868, CEP 01506-000, 
Sao Paulo, Brazil 

email: diego.goncalves@cruzeirodosul.edu.br 

^Astronomy Department, University of Wisconsin, Madison, 475 N. Charter St., WI 53711, USA 

^Department of Physics and Astronomy, the University of Western Ontario, London, Ontario, Canada, 
N6A 3K7 



- 2 - 



observed line profiles can provide the spectral index and the energy transfer rate 
of turbulence. Also, comparing two similar simulations with different viscous 
coefficients it was possible to constrain the turbulent cut-off scale. This may 
especially prove useful since it is believed that ambipolar diffusion may be one 
of the dominant dissipative mechanism in star-forming regions. In this case, the 
determination of the ambipolar diffusion scale may be used as a complementary 
method for the determination of magnetic field intensity in collapsing cores. We 
discuss the implications of our findings in terms of a new approach to magnetic 
field measurement proposed by Li & Houde (2008). 

Subject headings: ISM: magnetic fields, ISM: kinematics and dynamics, tech- 
niques: radial velocities, methods: numerical, statistical 

1. Introduction 

The interstellar medium (ISM) is known to be composed by a multi-phase, turbulent 

and magnetized gas (see Brunt & Heyer 2002, Elmegreen & Scalo 2004, Crutcher 2004, 
McKee & Ostriker 2007). However, the relative importance of turbulence and the magnetic 
field in the ISM dynamics and in the formation of structures is still a matter of debate. More 
specifically, typical molecular clouds present densities in the range of 10^ — 10^ cm~^, sizes 
L ~ 0.1 — lOOpc, temperature T ~ 10 — 20 K, and lifetimes that arc larger than the Jeans 
gravitational collapse timescale. The role of the magnetic field in preventing the collapse 
is hotly debated in the literature (see Pledge & Pudritz 2000, Palceta-Gongalves, de Juli & 
Jatenco-Pereira 2003, MacLow & Klessen 2004). Magnetic field can be removed from clouds 
in the presence of ambipolar diffusion arising from the differential drift of neutrals and ions 
(Mestel & Spitzer 1986, Shu 1983) and reconnection diffusion which arises from fast magnetic 
reconnection of turbulent magnetic field (Lazarian 2005). Nevertheless, the role of magnetic 
fields in the dynamics of ISM is difficult to underestimate. 

We feel that a lot of the unresolved issues in the theory of star formation are in part due 
to the fact that the amount of information on magnetic fields obtainable through presently 
used techniques is very limited. For molecular cloud the major ways of obtaining information 
about magnetic fields amount to Zeeman broadening of spectral lines, which provides mea- 
sures of the field strength along the line of sight (see Crutcher 1999) and the Chandrashekar- 
Permi (CP) method, which uses the statistics of polarization vectors to provide the amphtude 
of the plane of the sky component of the field (see Hildebrand 2000). However, since Zeeman 
measurements are restricted to rather strong magnetic fields (due to current observational 
sensitity) and therefore the measurements are restricted to dense clouds and new measure- 
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ments require a lot of observational time. At the same time, the CF method rehes on all 
grains being perfectly aligned, which is known not to be the case in molecular clouds (see 
Lazarian 2007 for a review). The CR technique is also known to systematically overestimate 
the field intensity (Houdc et al. 2009, Hildebrand et al. 2009), and to poorly map the mag- 
netic field topology for super- Alfvenic turbulence (see Falceta-Gongalves, Lazarian & Kowal 
2008). 

The difficulty of the traditional techniques call for new approaches in measuring as- 
trophysical magnetic fields. Recently, a number of such techniques has been proposed. For 
instance, Yan & Lazarian (2006, 2007, 2008) discussed using the radiative alignment of atoms 
and ions having fine or hyperfine split of the ground of metastable levels. The technique 
is based on the successful alignment of atoms in the laboratory conditions, but it requires 
environments where radiative pumping dominates the coUisional de-excitation of the levels. 

Another new approach which we dwell upon in this paper is based on the comparison 
of the ion-neutral spectral fines. Houde et al. (2000a, 2000b) identified the differences of 
the width of the lines of neutral atoms and ions as arising from their differential interaction 
with magnetic fields. It was assumed that because ions are forced into gyromagnetic motions 
about magnetic field lines that their spectral line profiles would thus reveal the imprint of 
the magnetic field on their dynamics. 

In particular, as observations of HCN and HCO+ in molecular clouds revealed signifi- 
cantly and systematically narrower ion lines, Houde et al. (2000a) proposed a simple expla- 
nation for these observations. The model was solely based on the strong Lorentz interaction 
between the ion and the magnetic field lines, but also required the presence of turbulent mo- 
tions in the gas. More precisely, it was found that the observations of the narrower HCO^ 
lines when compared to that from the coexistent HCN species could potentially be explained 
if neutral particles stream pass magnetic field lines with the entrained ions. Such a picture 
could be a particular manifestation of the ambipolar diffusion phenomenon. 

Although this model was successful in explaining the differences between the velocities 
of ions and neutrals, the quantitative description of the model of drift was oversimplified. 
For example, it was neither possible to infer anything about the strength of the magnetic 
field nor was the "amount" of ambipolar diffusion, which is at the root of the observable 
effect described by the model, quantifiable in any obvious manner. The main reason for 
these shortcomings resides in the way that turbulence and its interplay with the magnetic 
field were treated in the analysis of Houde et al. (2000a); a more complete and powerful 
model was required. The next step in the study of magnetized turbulence and ambipolar 
diffusion through the comparison of the coexistent ion/neutral spectral lines was taken by Li 
& Houde (2008) where a model of turbulence damping in partially ionized gas was employed. 
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Assumption that the main damping mechanism is associated with ambipolar diffusion, they 
deduced proposed a way for evaluating the strength of the plane-of-the-sky component of 
the magnetic field in molecular clouds. 

The main idea is that, in a magnetically dominated scenario, cloud collapse and magnetic 
energy removal may be accelerated due to ambipolar diffusion of ions and neutral particles. 
Although as gravity becomes dominant the collapsing cloud continuously drags material to 
its core including the ions, which are frozen to the magnetic field lines, magnetic pressure 
slows down their infall, but not that of the neutrals. At this stage most of the matter, in 
neutral phase, continues to decouple from the ionic fiuid and the field lines leading to the 
diffusion of magnetic energy and the collapse may develop further. This ion-neutral drift, 
excited by the ambipolar diffusion, is also responsible for damping the ion turbulent motions. 
The increase in the net viscosity of the flow provides a cut-off in the turbulent cells with 
turnover timescales lower than the period of collisions Ti,„ (see Lazarian, Vishniac & Cho 
[2004] for detailed review). 

Since the turbulent cascade is dramatically changed by the decoupling of the ion and 
neutral fluids the observed velocity dispersion could reveal much of the physics of collapse 
length scales. The interpretation and rehabihty of this technique, however, still need to 
be corroborated with more detailed theoretical analysis, as well as numerical simulations of 
magnetized turbulence. 

For the past decade, because of its complicated, fully non-linear and time-dependent na- 
ture, magnetized turbulence has been mostly studied by numerical simulations (see Ostriker, 
Stone &: Gammie 2001, Cho & Lazarian 2005, Kowal. Lazarian & Beresniak 2007). Sim- 
ulations can uniquely provide the three-dimensional structure for the density, velocity and 
magnetic fields, as well as two dimensional maps that can be compared to observations (e.g. 
column density, line profiles, polarization maps). Therefore, direct comparison of observed 
and synthetic maps may help reveal the magnetic topology and velocity structure. 

In this paper we re-examine the assumptions made in this model and test some of these 
assumptions using the MHD numerical simulations. In particular, we provide a number 
of numerical simulations of MHD turbulent flows, with different sonic and Alfvcnic Mach 
numbers. In §2, we describe the NIDR technique for the determination of the damping scales 
and magnetic fleld intensity from dispersion of velocities and the main theoretical aspects 
of turbulence in partially ionized gases. In §3, we describe the numerical simulations and 
present the results and the statistical analysis of the data. In §4, we discuss the systematic 
errors intrinsic to the procedures involved, followed by the discussion of the results and 
summary, in §5. 
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2. Turbulence in Partially Ionized Gas 



2.1. Challenge of interstellar turbulence 



In 1941, Kolmogorov proposed the well-known theory for energy cascade in incom- 
pressible fluids. Under Kolmogorov's approximation, turbulence evolves from the largest to 
smaller scales, up to the dissipation scales, as follows. Within the so-called inertial range, 
i.e the range of scales large enough for dissipation to be negligible but still smaller than the 
injection scales, the energy spectrum may be well described by. 



where e is the energy transfer rate between scales and a ~ 5/3. In this approximation, 
within the inertial range, the energy transfer rate is assumed to be constant for all scales. 
Therefore, integrating Eq. 3 over /c for a = 5/3, we obtain. 



However, reality is far more complicated. First, the ISM is threaded by magnetic flelds, 
which may be strong enough to play a role on the dynamics of eddies and change the 
scaling relations. Second, observations suggest that the ISM is, at large scales, is highly 
compressible. Third, many phases of the ISM (see Draine & Lazarian 1999 for typical 
parameters) are partially ionized. 

Attempts to include magnetic flelds in the picture of turbulence include works by Irosh- 
nikov (1964) and Kraichnam (1965), which were done assuming that magnetized turbulence 
stays isotropic. Later studies proved that magnetic fleld introduces anisotropy into turbu- 
lence (Shebalin, Matthaeus, & Montgomery 1983, Higdon 1984, Zank & Matthaeus 1992, 
see also book by Biskamp 2003). 

Goldreich & Sridhar (1995, henceforth GS95) proposed a model for magnetic incom- 
pressible turbulenc^ based on the anisotropies in scaling relations, as eddies would evolve 
differently in directions parallel and perpendicular to the fleld lines: V^/Ay ~ v\/\, where 



^In the original treatment of GS95 the description of turbulence is limited to a situation of the velocity 
of injection at the injection scale Vl being equal to Va- The generalization of the scalings when Vl < Va 
can be found in Lazarian & Vishniac (1999). The generalization for the Vl > Va is also straightforward (see 
Lazarian 2006). 



(1) 




(2) 
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A|j is the parallel scale of the eddy and A is its perpendicular scale. These scales are measured 
in respect to the /oca^ magnetic field. Combining this to the assumption of self-similarity 
in energy transfer rate, we get a Kolmogorov-like spectrum for perpendicular motions with 
a = 5/3 and, most importantly, the anisotropy in the eddies scales as Ay oc A^/^. 

In spite of the intensive recent work on the incompressible turbulence (see Boldyrev 
2005, 2006, Beresnyak & Lazarian 2006, 2009), we feel that the GS95 is the model that can 
guide us in the research in the absence of a better alternative. The generalization of the 
GS95 for compressible motions are available (Lithwick & Goldreich 2001, Cho & Lazarian 
2002, 2003) and they consider scalings of the fast and slow MHD modes. 



2.2. Turbulence damping in partially ionized gas 

As stated before, the turbulent cascade is expected to develop down to 
scales where dissipation processes become dominant. The dissipation scales are 
associated to the viscous damping, which is responsible for the transfer of kinetic 
into thermal energy of any eddy smaller than the viscous cutoff scale. 

In the ISM, e.g. in cold clouds, the gas is partially ionized and the coupling 
between neutrals, ions and magnetic fields gives rise to interesting processes. As 
far as damping is concerned one of the most interesting is the energy dissipation 
as the motions of ions and neutral particles decouple. While the issue of turbu- 
lence dissipation has been discussed extensively in the literature (see Minter & 
Spangler 1997), a generalization of the GS95 model of turbulence for the case 
of the partially ionized gas was presented in Lazarian, Vishniac &z Cho (2004) 
(LVC04). In their model, if the eddy turnover time (r) gets of order of the 
ion-neutral collision rate (t^^) two fluids are strongly coupled. In this situation 
a cascade cut-off is present. 

In a strongly coupled fluid, using the scaling relation for the inertial range 
^^damp ~ f/inj(^in}^damp)^/^ the damping scale is given by (LVC04), 

^damp ~ Amfp ( J ( ] /„, (3) 



^The latter issue does not formally allow to describe turbulence in the Fourier space, as the latter calls 
for the description in respect to the global magnetic field. 
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where /„ is the neutral fraction, Amfp and the mean free path and sound speed 
for the neutrals, respectively, Va is the Alfven speed, and subscript "inj" refers 
to the injection scale. Since the Alfvn speed depends on the magnetic field, 

Va = -B(47rp)"^/2, Eq.(3) is rewritten as: 

Therefore, for a given molecular cloud, if the decoupling of ions and neutrals 
is the main process responsible for the ion turbulence damping Eq.(4) may be 
a complementary estimation of B. The main advantage of this method is that 
B is the total magnetic field and not a component, parallel or perpendicular to 
the LOS, as respectively obtained from Zeeman or CF-method from polarization 
maps. 



2.3. Approach by Li &; Houde 2008 

From the perspective of the turbulence above we can discuss the model adopted by Li & 
Houde (2008) for their study. The authors considered that damping of ion motions happen 
earlier than those by neutrals at sufficiently small scales. At large scales, ions and neutrals 
are well coupled through flux freezing and their power spectra should be similar. At small 
scales the ion turbulence damps while the turbulence of neutral particles continues cascading 
to smaller scales. This difference may be detected in the velocity dispersions (a) obtained 
from the integration of the velocity power spectrum over the wavenumber k, 



a^k) oc / P{k')dk' 

Jk 

~ bk-"", (5) 

considering P{k) a power-law spectrum function. Since the turbulence of ions is damped at 
the diffusion/dissipation scale (^d), while the turbulence of neutral particles may develop 
up to higher wavenumbers we may consider that the ions and neutral particles present the 
same distribution of velocities (well coupled) for L > Lj^. In this sense, the dispersion of 
neutral particles may be written as. 



/■oo 

al{k) oc / P{k')dk' 

Jk 
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rkj3 POO 

= / P{k')dk'+ / P{k')dk' 

Jk Jko 

~ aUk) + bk^^, (6) 

where A;d ~ L^^. Eq.(5) may be directly compared to Eq. 2. In this case, we would obtain 
n ~ 2/3 and b is related to the energy transfer rate e. Therefore, once the fitting parameters 
of Eq. 1 are obtained from the observational data, it is possible to obtain the cascading 
constants e and a. 

With the dispersion of velocities for both ions and neutrals at different scales k it is 
possible to calculate the damping scale /cd- Prom Eq.(5), the dispersion of neutral particles 
provides b and n constants and, by combining neutral and ion dispersions, it is possible to 
get /cd (Eq.6), i.e. the damping scale. 

Finally, as proposed by Li & Houdc (2008) in a different context, it is possible to to 
evaluate magnetic field strength by Eq.(4). We feel that the procedure of magnetic field study 
requires a separate discussion, due to its complexity, but in what follows we concentrate on 
the interesting facts of observational determining of the characteristics of turbulence and its 
damping that are employed in the technique by Li & Houde (2008). 



2.4. Observational perspective 

Eqs. (1) and (2) are based on the dispersion of a three-dimensional velocity field, i.e. 
subvolumes with dimensions k"^. Observational maps of line profiles, on the other hand, 
provide measurements of the velocity field integrated along the fine of sight (LOS) within the 
area of the beam, i.e. a total volume of k'^k^l^ (with k = l/l and /cmin = ^/L, as L represents 
the total depth of the structure observed - typically larger than /). Also, velocity dispersions 
are obtained from spectral line profiles, which are strongly dependent on the column density, 
i.e. the distribution of matter along the LOS. These factors make the comparison between 
observed lines and theoretical distribution of velocity fields a hard task. 

Fortunately, 3-dimensional numerical simulations of MHD turbulence may be useful in 
providing both the volumetric properties of the plasma parameters as well as their synthetic 
measurements projected along given lines of sight, which may be compared directly to ob- 
servations, such as the spectral line dispersion. In this sense, based on the simulations of 
Ostriker et al. (2001), Li & Houde (2008) stated that the actual dispersion of velocity is, 
approximately, the minimum value of the LOS dispersion, at each beam size I, obtained in a 
large sample of measurements. However, Ostriker et al. (2001) presented a single simulation, 
exclusively for supersonic and sub-alfvenic turbulent regime, with limited resolution (256^). 
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They also did not study increased viscosity, nor the correlation of minima of the synthetic 
dispersion and the turbulent regimes and the distribution of gas along the LOS. 

In the following sections we will describe the details regarding the estimation of line 
dispersions, but now comparing it with a larger set of numerical simulations with different 
turbulent regimes and with finer numerical resolution. The idea is to determine whether 
the technique is useful or not, and if there is any limitations with the different turbulent 
regimes. These tests are mandatory to ensure the applicability of the NIDR method to ISM 
observations. 



In order to test the NIDR model, i.e to verify if the minimum dispersion of the velocity 
measured along the line of sight for a given beamsize / x / is aproximately the actual value 
calculated for a volume we used a total of 12 3-D MHD numerical simulations, with 512'^ 
resolution, for 6 different turbulent regimes as described in Table 1, but repeated for viscous 
and inviscid models. 

The simulations were performed solving the set of ideal MHD isothermal equations, in 
conservative form, as follows: 



3. 



Numerical Simulations 



dp 
dt 



+ V • (pv) = 0, 



(7) 



dpv 



dt 



+ V 



pvv + {p + 




BB 



f, 



(8) 



dB 
'dt 



V X (v X B) = 0, 



(9) 



V-B = 0, 



(10) 



P = c^p, 



(11) 



where p, v and p are the plasma density, velocity and pressure, respectively, B = V x A 
is the magnetic field, A is the vector potential and f = fturb + fvisc represents the external 
source terms, responsible for the turbulence injection and explicit viscosity. The code 
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solves the set of MHD equations using a Godunov-type scheme, based on a 
second-order-accurate and the non-oscillatory spatial reconstruction (see Del 
Zanna et al. 2003). The shock-capture method is based on the Harten-Lax-van 
Leer (1983) Riemann solver. The magnetic divergence-free is assured by the 
use of a constrained transport method for the induction equation and the non- 
centered positioning of the magnetic field variables (see Londrillo &: Del Zanna 
2000). The code has been extensively tested and successfully used in several 
works (Falceta-Gongalves, Leizcirian &: Kowal 2008; Leao et al. 2009; Burkhcirt 
et al. 2009; Kowal et al. 2009; Falceta-Gongalves et al. 2010). 

The turbulence is triggered by the injection of solenoidal perturbations in Fourier space 
of the velocity field. Here, we solve the explicit viscous term as fvisc = — pi^V^v, where u 
represents the viscous coefficient and is set arbitrarily to simulate the increased viscosity of 
the ionic flows due to the ambipolar diffusion. We run all the initial conditions given in Table 
1 for both u = and u = 10~^, representing the neutral and ion particles fluids, respectively. 
Each simulations is initiated with an uniform density distribution, threaded by an uniform 
magnetic field. The simulations were run until the power spectrum is fully developed. The 
simulated box boundaries were set as periodic. 

In Fig. 1 we show the resulting velocity power spectra of four of our models, representing 
the four different turbulent regimes, i.e. (sub)supersonic and (sub)super-Alfvcnic. Solid 
lines represent the non-viscous cases, and the dotted line the viscous cases. The spectra are 
normalized by a Kolmogorov power function P oc A;^^/"^. The incrtial range of the scales 
is given by the horizontal part of the spectra. For the inviscid fluid, subsonic turbulence 
presents approximately flat spectra for 2 < A; < 50. Supersonic turbulence, on the other 
hand, shows steeper power spectra within this range. Actually, as shown from numerical 
simulations by Kritsuk et al. (2007) and Kowal & Lazarian (2007), shocks in supersonic 
flows are responsible for the fllamentation of structures and the increase in the energy flux 
cascade, resulting in a power spectrum slope ~ —2.0. For k > 50, the power spectra show 
a strong damping of the turbulence, resulting from the numerical viscosity. For the viscous 
fluid, the damped region is broadened (/Ccutoff ~ 20), due to the stronger viscosity. 

4. Relationship between 2D and 3D dispersion of velocities 

4.1. Compciring the synthetic to the 3-dimensional dispersion of velocities 

Theoretically, as given by Eq.(6), the difference between the two spectra for each run 
may be obtained from the observed dispersion of velocities. The next step then is to obtain 
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Fig. 1. — Velocity power spectra of four of the models described in Table 1, separated by 
turbulent regime. The spectra are normalized with a Kolmogorov power function (P oc 
k~^^^). Solid lines represent the non- viscous fluid, and the dotted line the viscous fluid. 



the dispersion of velocity, for different scales /, from our simulations. However, as explained 
previously, there are two different methods to obtain this parameter. One represents the 
actual dispersion, calculated within subvolumes of the computational box, while the second 
represents the observational measurements and is the dispersion of the velocity within the 
subvolume I'^L (assuming the gas is optically thin), where L is the total depth of the box. 
In order to match our calculations to observational measurements we will use the density 
weighted velocity v* — pv (see Esquivel & Lazarian [2005]), which characterizes the line 
emission intensity proportional to the local density. 

To obtain the actual dispersion of velocities as a function of the scale I, we subdivide the 
box in N volumes of size l^. Then, we calculate the dispersion of v*, normalized by the sound 
speed Cs, as the mean value of the local dispersions obtained for each subvolume. For the 
synthetic observational dispersion, we must firstly choose a given line of sight (LOS). Here, for 
the results shown in Fig. 2 we adopted x-direction. After, we subdivide the orthogonal plane 
(y-z) in squares of area representing the beamsize. Finally, we calculate the dispersion of 
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V* within each of the volumes PL, normahzed by the sound speed c^, for different values of 
l/L. The results of these calculations for each of the non- viscous models, is shown is Fig. 2. 

The solid line and the triangles represent the average of the actual mean dispersion 
of velocities, while the crosses represent each of the synthetic observed dispersion within 
PL. Regarding the synthetic observational measurements, we see that increasing I results 
in a decrease in the dispersion, i.e. range of values, of (7„*. Also, since we use the density 
weighted velocity v*, the mass distribution plays an important role in the calulation of 
a"„*. Denser regions will give a higher weight for their own local velocities and, therefore, if 
several uncorrelated denser regions are intercepted by the LOS, ay. will probably be larger. 
Therefore, we may understand the minimum value of a^* (/) as the dispersion obtained for 
the given LOS that intercepts the lowest number of turbulent sub-structures. If a single 
turbulent structure could be observed, then a^* (l) would tend to the actual volumetric value 
if the overdense structure depth is ~ Z. Also, as you increase I, the number of different 
structures intercepting the line of sight increases, leading to larger values of the minimum 
observed dispersion. On the other hand, the maximum observed dispersion is directly related 
to the LOS that intercepts most of the different turbulent structures. Since this number 
is unlikely to change, the maximum observed dispersion decreases with / simply because 
of the larger number of points for statistics. As / — )■ L, ay*{l) gets closer to the actual 
volumetric dispersion. However, as noted in Fig. 2, the obtained values for I — )■ L 
are slightly different. This is caused by the anisotropy in the velocity field 
regairding the magnetic field, as the velocity components may be different along 
and perpendiculeir to B. 

Despite of this effect, the results presented in this work do not change when a 
different orientation for the line of sight is chosen. Even though not shown in Fig. 
2, we have calculated the dispersion of velocity for LOS in y and z-directions. 
The general trends shown in Fig. 2 ctre also observed, but a slight difference 
appecirs as / — )■ L, exactly as explained above. This difference is expected to be 
seen in sub-alfvenic cases because of the anisotropy in the velocity distribution. 

It is clear from Fig. 2 that the actual dispersion of velocities and a given observational 
line- width may be very different. Li & Houde (2008), based on Ostriker et al.'s work, assumed 
that if one chooses, from a large number of observational measurements along different LOS's, 
the minimum observational dispersion as the best estimation for the actual dispersion, the 
associated error is minimized. Considering the broad range of observed dispersions obtained 
from the simulations for a given /, the minimum value should correspond to the actual 
velocity dispersion. Actually, from Fig. 2 we see that the vahdity of such statement depends 
on I and on the turbulent regime, though as a general result the scaling of the minimum 
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observed dispersion follows the actual one. 

For the subsonic models, the actual dispersion is lower than cr^j(/), with increasing 
difference as Z/L — > 0. In these models, we see that for I > 0.05L there is a convergence 
of the actual dispersion to the synthetic observational measurements. At these scales the 
minimum value of crv*{l) is a good estimate of the velocity dispersion of the turbulence at 
the given scale I. The difference between both values is less than a factor of 3 for all Vs, 
being of a few percent for / > 0.03L. In this turbulent regime, mainly at the smaller scales, 
(Tv*{l) overestimates the true dispersion. For larger scales the associated error is very small 
and the two quantities give similar values. 

On the other hand, for the highly supersonic models (M5 ~ 7.0), the minimum value of 
a„j(/) underestimates the actual value, at most scales. Under this regime, the difference to 
the actual dispersion is of a factor ~ 2—4. The best matching between the two measurements 
occured for the marginally supersonic cases (M5 ~ 1.5). 

As a major result we found that the uncertainties associated to the NIDR technique 
depend on the sonic Mach number of the system, though the associated errors are not extreme 
in any case. We see no major role of the Alfvenic Mach number on this technique. 

4.2. Minimum 2D velocity dispersion vs 3D statistics 

We have shown that the synthetic observed dispersion minima represent a fair approxi- 
mation for the actual 3-dimcnsional dispersion of velocities for the supersonic models, though 
it is slightly overestimated in subsonic cases, and underestimated in highly supersonic cases. 
What is the physical reason for that? 

One of the most dramatic differences between subsonic and supersonic turbulence is 
the mass density distribution. Subsonic turbulence is almost incompressible, which means 
that density fluctuations and contrast are small. Supersonic turbulence, on the other hand, 
present strong contrast and large fluctuation of density within the volume. Strong shocks 
play a major role on the formation of high density contrasts, and is the main cause of high 
density sheets and filamentary structures in simulations. Strong shocks also modify the 
turbulent energy cascade, opening the possibility for a more efficient transfer of energy from 
large to small scales, resulting in a steeper energy spectra (typically with index ~ —2, instead 
of the Kolmogorov's ~ —5/3). 

Observationally, the determination of the velocity dispersion along the hue of sight is 
always biased by the density distribution, i.e. a given line profile depends on both the emis- 
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sion intensity and the Doppler shifts. For an optically thin gas, the emission intensity is 
directly dependent on the density of the gas. Compared to subsonic, we expect the super- 
sonic turbulence to present a larger contrast of density in structures. Because of that, the 
sizes and the number of structures intercepting the line of sight may have a deep impact in 
the determination of the observational dispersion of lines. In order to check this hypothesis 
we calculated the number of structures, and their average sizes, along each of the synthetic 
lines of sight of the cube and determined the correlation with the velocity dispersion. The 
number and depths of the structures were obtained by an algorithm that identi- 
fies peeiks in density distributions. Bcisically, the algorithm follows three steps. 
Firstly, for each line of sight with beamsize Z^, it identifies the maximum pecik 
of density and uses a threshold defined as the half value of this maximum of 
density. Secondly, it removes all cells with densities lower than the selected 
threshold. The remaining data represents the dominant structures within the 
given line of sight. Finally, the algorithm follows each line of sight detecting 
the discontinuities, created by the use of a threshold, and calculates the sizes of 
these structures. 

In Fig. 3 we show the correlation between the synthetic velocity dispersion and the 
number of structures intercepted in all LOS, in x-direction, for our cube of the Model 3. 
For all models the result is very similar, though not shown in these plots. As noticed, the 
minimum dispersion corresponds to the LOS in which there is only one structure intercepted, 
i.e. there is only one source that dominates the emission line. It makes complete sense if this 
single source is small in depth and, as a volume, we have the dispersion of a volume ~ l^. As 
we increase the number of structures intercepted by the LOS, each one contributes with a 
different Doppler shift, resulting in a larger dispersion. However, as we can see from Fig. 3, 
the maximum dispersion also corresponds to a single structure in the LOS. The reason is that 
we calculate the threshold for capturing clumps with the FWHM of the highest peak of the 
given LOS. If the LOS intercepts "voids", which are typically very large compared to clumps, 
the algorithm results in a number of structures equal one but its depth, and consequently 
also its velocity dispersion, is large. Taking into account observational sensitivity, these low 
density regions are irrelevant. Furthermore, this picture is also useful to understand the 
scaling relation of As we increase I, increasing the beamsize, the number of structures 
in the LOS is higher resulting in the increase in velocity dispersion. 

In Fig. 4, we compare the average number of structures intercepting the LOS's and their 
average sizes as a function of the sonic Mach number. It is shown that both the number 
of sources and their intrinsic sizes are inverselly correlated with the sonic Mach number. 
Subsonic models present lower contrast of densities, which corresponds to larger overdense 
structures. For Ms = 0.7 we obtained a range of densities 0.3 < p/ p < 3. As discussed 
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above, large structures correspond to integrated volumes ~ PL, which deviates from the 
actual 3D velocity dispersion. The result is that the minimum synthetic dispersion obtained 
from subsonic turbulence will overestimate the actual value. For supersonic models, where 
clumps are systematically small, the dispersion minima correspond, from Fig. 3, to the 
single structures with volumes ~ very close to the actual values. On the other hand, for 
/ L, the result is underestimated. Here, for Ms = 1.5 we obtained a range of densities 
0.08 < p/p < 10, while for Mg = 7.0 we obtained 0.01 < p/p < 90. 



4.3. Turbulence dissipation scales 

In order to obtain the dissipation scales of the ionized flows, we must apply Eqs. (5) 
and (6) and 2 to the simulated data. For the inviscid flows, we calculate b and n (Eq.5) 
using both the actual and minimum observed velocity dispersions. Then, we use the data 
from the viscous simulation to calculate the difference a^ — af. Finally, from Eq.(6), knowing 
(7^ — crf, n and h it is possible to obtain k-Q- In Fig. 5 we show the data used to compute 
Eqs. (5) and (6), i.e. both the actual (lines) and the minimum observed (symbols) dispersions 
of velocity for the viscous (stars) and inviscid fluids (squares). The fit parameters (Eq.5) 
for the inviscid simulated data arc shown in Table 2, where bac and n^c where obtained for 
the actual velocity dispersion and bobs and Uobs for the synthetic observed line widths. We 
see that n increases with the sonic Mach number (M5), as explained below. In Table 2, we 
also present the damping scale Ld, obtained from Eq.(6). In the last column we show the 
ratio between the actual and observational scales L^/L'^^^. The ratio of the obtained scales 
showed a maximum difference of a factor of 5 between the actual value of the dispersion and 
the one obtained from the synthetic observational maps. Also, compared to the expected 
values obtained visually from the spectra (Fig.l), there is a good correspondence with L°^^ 
given in Table 2. This fact shows that the method indeed might be useful. 

Regarding the spectral index a (Eq.l), Table 2 gives a ~ 1.4 — 1.8 for the subsonic 
models, and a ~ 2.0 — 2.1 for the supersonic models. These parameters may be directly 
compared to the values a ~ 1.7 and a ~ 2.0 expected for theoretical incompressible and 

compressible turbulent spectra, respectively. Also, the increase of h as Ms increases shows 
that the energy transfer rate is larger for compressible models. From Table 2, we obtain an 
averaged value of e ~ 0.6 for subsonic and 1.4 for supersonic models, in code units. 



- 16 - 



Table 1: Description of the simulations - .B is assumed in x-direction 



Model Ms Ma Description 



1 


0.7 


0.7 


subsonic & sub-Alfvenic 


2 


1.5 


0.7 


supersonic & sub-Alfvenic 


3 


7.0 


0.7 


supersonic & sub-Alfvenic 


4 


0.7 


7.0 


subsonic & super-Alfvenic 


5 


1.5 


7.0 


supersonic & super-Alfvenic 


6 


7.0 


7.0 


supersonic & super-Alfvenic 



sonic Mach number {Ms = {v/cs)) 
Alfvenic Mach number (M^ = (v/Va)) 



Table 2: Parameters of best fit and damping scales 



Ms 


Ma 


bac 


actual 

nac 


7" ac 


bobs 


synthetic maps 

Hobs 




r ac / r obs 


0.7 


0.7 


0.63(16) 


0.39(2) 


0.244(21) 


1.01(9) 


0.80(5) 


0.088(13) 


2.7 


0.7 


7.0 


1.38(12) 


0.65(3) 


0.246(21) 


1.00(10) 


0.76(4) 


0.215(17) 


1.1 


7.0 


0.7 


1.58(21) 


1.01(4) 


0.162(17) 


1.58(13) 


1.08(5) 


0.318(23) 


0.5 


7.0 


7.0 


1.82(18) 


1.05(4) 


0.041(6) 


2.51(18) 


0.95(6) 


0.197(25) 


0.2 



and L°^^ are placed in terms of the total size of the box, i.e. Ld/L. 
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4.4. Accuracy of the method 

Before describing a potential technique for future observations, we must stress out that 
this work may be divided in two independent parts. The first is related to the characterization 
of the energy power spectrum of turbulence, including the energy transfer rate and cut-off 
length, while the second is based on the use of this damping length which is used in the Li 
& Houde (2008) approach to determine the magnetic field. 

As mentioned before, in magnetized partially ionized gases, several different mechanisms 
may act in order to damp turbulent motions. The determination of the magnetic field 
intensity from the damping scale is strictly dependent on the assumption that the turbulent 
cut-off is due to the ambipolar diffusion, i.e. the ion-neutral diffusion scale is larger than 
the scales of any other dissipation mechanism. Unfortunately, this could not be tested in 
the simulations since we did not include explicit two fluid equations to check the role of 
ambipolar diffusion and other damping mechanisms in the turbulent spectra. 

In the simulations, the cut-off is obtained via an explicit and arbitrary viscous coefficient. 
The result is clear in the power spectra (Fig. 1), where the damping length shifts to lower 
values of k. Prom those, we found out that the "observational" velocity dispersions for a 
given beam size Ixl, in most cases, do not coincide with the actual dispersion at scale /. Also, 
the estimation from the minimum value of the observed dispersions may also be different 
from the expected measure. The associated errors depend on the sonic Mach number and 
on the scale itself, as shown in the previous section. However, the parameters obtained for 
the power-law of the synthetic observations and actual 3-dimensional distributions showed 
to be quite similar. The theory behind this method is very simple, but still reasonable, as it 
assumes that the two fiuids (ions and neutrals) would present the same cascade down to the 
dissipation scale, when they decouple, where the ion turbulence would be sharply damped. 
We believe that these errors may, eventually, be originated by the short inertial range of the 
simulated data. We see from Fig. 1 that the turbulent damping range is broad, and not 
sharp compared to the inertial range, as assumed in this model. In the real ISM, the power 
spectrum presents a constant slope within a much broader range, and the errors with real 
data may be smaller. 



5. Discussion 

In this work we studied the relationship between the actual 3-D dispersion of velocities 
and the ones obtained from synthetic observational line profiles, i.e. the density weighted 
line profile widths, along different LOS's. We study the possibility of the scaling relation of 
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the turbulent velocity dispersion being determined from spectral line profiles. If correct, the 
observed line profiles could allow us to determine in details the turbulence cascade and the 
dissipation lengths of turbulent eddies in the ISM. However, in order to check the validity 
of this approach, we performed a number of higher resolution turbulent MHD simulations 
under different turbulent regimes, i.e. for different sonic and Alfvenic Mach numbers, and 
for different viscosity coefficients. 

Based on the simulations, we showed that the synthetic observed line width (a^) is 
related to the number of dense structures intercepted by the LOS. Therefore, the actual 
dispersion at scale / tends to be similar to the line width obtained by a LOS within a 
beam size / x / intercepting a single dense structure. It means that, the minimum observed 
dispersion may be the best estimative for the actual dispersion of velocities, if a large number 
of LOS is considered, as assumed in the theoretical model. 

Moreover, by adopting a power-law for the spectrum function of ay, we were able to 
estimate the spectral index n and constant b (Eq. 5), given in Table 2. We see a good 
correspondence between the parameters obtained for the 3-D dispersions and synthetic line 
profiles. Furthermore, since n is associated with the turbulent spectral index a, and b with 
the energy tranfer rate between scales, line profiles may be useful in characterizing the 
turbulent cascade. 

Also, we showed that the models under similar turbulent regime but with different 
viscosities will result in different dispersions of velocity, on both 3-D and synthetic line 
profile measurements. The difference of ay for the inviscid and viscous models, associated 
with the parameters n and b previously obtained, gives an estimative of the damping scale 
Ld of the viscous model (Eq. 6). The ratio between the damping scales obtained from the 
3-D and synthetic profile dispersions vary only by a factor < 5. 

A good estimate of the cut-off scale of the ISM turbulence may bring light to much of 
the uncertainties about the mechanisms that are responsible for the damping of the turbulent 
eddies. As we showed, if ambipolar diffusion is the main phenomenon responsible for the 
dissipation of turbulence, then it is possible to provide another method for the determination 
of magnetic field strength in dense cores, besides Zeeman splitting. 

We believe that, i - the numerical resolution used in our models, and ii - the single 
fluid aproximation, with different viscosities to simulate neutrals and ions separately, are the 
main limitations of this work. The ISM may present more than 5 decades of inertial range 
in its power spectrum while numerical simulations are, at best, limited to 2-3. On the other 
hand, since the validity of the theoretical approximation presented in this work depends on 
the broadness of the inertial range, we expect this model to be even more accurate for finer 
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resolution, or for the ISM itself. Under the single fluid aproximation made in this work we 
fail in correlating the ambipolar diffusion with increased viscosity and, therefore it is not 
possible to test the estimation of B from the damping scales as proposed above, though 
the results related to the damping length and energy transfer rates remain unchanged. In 
a two-fluid simulation this is more likely to be fulfilled. We are currently implementing the 
two-fiuid set of equations in the code, and will test this hypothesis in the future. 

6. Summary 

In this work we presented an extensive analysis of the applicability of the NIDR method 
for the determination of the turbulence damping scales and the magnetic field intensity, if 
ambipolar diffusion is present, based on numerical simulations of viscous MHD turbulence. 
We performed simulations with different characteristic sonic and Alfvenic Mach numbers, 
and different explicit viscous coefficients to account for the physical damping mechanisms. 
As main results we showed that: 

• the correspondence between the synthetic observational dispersion of velocities (i.e. 
from the 2D oscrvcd maps) and the actual 3-dimensional dispersion of velocities de- 
pends on the turbulent regime; 

• for subsonic turbulence, the minimum inferred dispersion tends to overestimate the 
actual dispersion of velocities for small scales {I « L), but presented good convergence 
at large scales (Z — >■ L); 

• for supersonic turbulence, on the other hand, there is a convergence at small scales 
{I « L), but the minimum inferred dispersion tends to underestimate the actual 
dispersion of velocities at large scales (Z — > L); 

• even though not precisely matching, the actual velocity and the minimum velocity 
dispersion from spectral lines were well fitted by a power-law distribution. We obtained 
similar slopes and linear coefficients for both measurements, with a ~ —1.7 and —2.0 
for subsonic and supersonic cases, respectivelly, as expected theoretically; 

• the damping scales obtained from the fit for the both cases were similar. The difference 
between the scales obtained from the two fits was less than a factor of 5 for all models, 
indicating that the method may be robust and used for observational data; 

The work presented in this paper tests some of the key assumptions important process 
by technique of Li & Houde (2008). Evidently, more work is still required in order to test the 
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fuU range of applicability of the method (e.g., test cases of both magnetically and neutral 
driven turbulence). The aforementioned implementation of two-fluid numerical simulations 
to better mimic ambipolar diffusion will be an important step in that direction. 
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Dispersion of Velocity versus Beam Size - X 
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Fig. 2. — Dispersion of velocity for models 1 — 6 obtained for different beam sizes (P) 
integrated along x-direction. The solid line (and triangles) represent the actual dispersion, 
defined as the averaged dispersion measured within all cubes with size For a given 
beamsize, the crosses represent each of the synthetic observed dispersion within l^L. 
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Fig. 3. — Correlation of the synthetic velocity dispersion and the number of structures 
intercepted in the line of sight. Each triangle correspond to all LOS's calculated from a cube 
of model 3. 
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Fig. 4. — Average number of structures within all LOS's computed for all models as a 
function of the sonic Mach number. 
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Fig. 5. — Actual (lines) and the minimum observed (symbols) dispersions of velocity for the 
viscous (stars) and inviscid fluids (squares). 



